Rare decay 7T° — ► e + e : theory confronts KTeV data. 
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Within the dispersive approach to the amplitude of the rare decay ir° — » e + e~ the nontrivial 
dynamics is contained only in the subtraction constant. We express this constant, in the leading 
order in (m e /A) 2 perturbative series, in terms of the inverse moment of the pion transition form 
factor given in symmetric kinematics. By using the CELLO and CLEO data on the pion transition 
form factor given in asymmetric kinematics the lower bound of the decay branching ratio is found. 
The restrictions following from QCD allow us to make a quantitative prediction for the branching 
B (tt° — > e + e~) = (6.2 ± 0.1) ■ 10 -8 which is 3cr below the recent KTeV measurement. We confirm 
our prediction by using the quark models and phenomenological approaches based on the vector 
meson dominance. The decays n —> l + l~ are also discussed. 
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1. Experimental measurements of neutral pseudoscalar meson decays into lepton pairs and the comparison with 
theoretical predictions offer an interesting way to study long-distance dynamics in the Standard Model. Recently, the 
KTeV E799-II experiment at Fermilab has observed ir° — » e + e~ events using Kt, — > 37r decay as a source of tagged 
neutral pions [l[ . The branching ratio of the pion decay into an electron-positron pair was determined to be equal to 
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where xo = (m e + e - /m^ ) 2 is the Dalitz variable. By extrapolating the full radiative tail beyond xd > 0.95 and 
scaling the result back up by the overall radiative correction 0, Q to find the lowest-order rate for ir° — > e + e~ the 
KTeV Collaboration obtained 
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^no-rad \~ 



(7.49 ±0.29 ±0.25) • 10" 
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The rare decay 7r° — > e + e~ has been studied theoretically over the years, starting with the first prediction of 
the rate by Drell Since no spinless current coupling of quarks to leptons exists, the decay is described in the 
lowest order of QED as a one-loop process via the two-photon intermediate state, as shown in Fig. 1. A factor of 
2 (m e /mn) corresponding to the approximate helicity conservation of the interaction and two orders of a suppress 
the decay with respect to the 7r° — > 77 decay, leading to an expected branching ratio of about 10 -7 . In the Standard 
Model contributions from the weak interaction to this process are many orders of magnitude smaller and can be 
neglected. The interaction of leptons and quarks with leptoquarks is a possible mechanism for the pion decay from 
physics beyond the Standard Model. The confrontation of the theory and experiment will have some influence on the 
problem of strong sector contribution to the muon anomalous magnetic moment g — 2 0, Q . 
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FIG. 1: Triangle diagram for it — > e + e process with a pion ir° — > 7*7* form factor in the vertex. 
To lowest order in QED the normalized branching ratio is given by 
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0.988 7]. In this article, we describe the process n° — > e + e by the 
diagram shown in Fig.l, where F! n . 7 » 7 . is the form factor of the transition 7r° — > 7*7* with off-shell photons. The 
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reduced amplitude A can be written as 



A (, 2 ) = I f d 4 7 q2 2 k2 \ ( f 5 r^rr (-* 2 , - (* - «) 2 ) - (4) 

rJ ^ (k 2 + ie) Uk - qf + ie\ Uk ~ pf - ml + ie\ V ' 

where q 2 = m 2 ,p 2 = m 2 . We put the sign minus in the arguments of the form factor explicitly to emphasize that 
Eq. (|4]) is written in the Minkowski space. The form factor is normalized as F 7r7 » 7 »(0, 0) = 1 and falls down quite 
rapidly in the Euclidean region of momenta to provide the ultraviolet convergence of the integral. A number of model 
calculations of the amplitude A (q 2 ) was performed 0, H, M, E3, EH EH, Eil by employing different shapes of the form 
factor F,r 7 . 7 *. We discuss some of them below. 

The aim of the present paper is to calculate the branching ratio B — > e + e _ ) and estimate the uncertainties by 
using the available experimental and theoretical information on the pion transition form factor. In particular, the 
important constraints follow from the results obtained by the CELLO and CLEO collaborations and restrictions set 
by QCD. 

2. First, we derive a suitable representation for the amplitude in Eq. ([4]) which would help us to perform a 
straightforward analysis by using the available information on the pion transition form factor. To do this, we employ 
the dispersive approach to the calculation of the amplitude developed in many papers (see, e.g. [12| and references 
therein) . The imaginary part of the amplitude in Eq. JH) 



comes from the contribution of real photons in the intermediate state and is model independent since i ? W7 » 7 * (0, 0) = 1. 
Using _4| 2 > (Im^4) 2 and neglecting radiative corrections one can get the well-known unitary bound for the branching 
ratio in Eq. © 

B (tt° -> e+e~) > B unitary (vr° -» e+er) = 4.69 • 1(T 8 . (6) 
A once-subtracted dispersion relation for the amplitude in Eq. ([4]) is written as [l2j 



A(^)=A( q 2 = 0) + q - r d s^4 y (7) 
7T Jo s i s ~ Q ) 



The second term in Eq. |7]) takes into account strong q 2 dependence of the amplitude around the point q 2 = 
occurring due to the branch cut coming from the two-photon intermediate state. Integrating Eq. ([7]) one arrives for 
q 2 > 4m 2 e at [13, LH, EBl 



ReA(q 2 ) =A(q 2 = 0) 



Pe(q 2 
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where Li2 (z) — — J* (dt/t)\n(l — t) is the dilogarithm function 1 . For the pion in the leading order in (m e /ra^) 2 , 
one gets 



ReA(ml) =A(q 2 = 0)+ In 2 f ^ j + ^. (9) 

Thus, the nontrivial dynamics is only contained in the subtraction constant A (q 2 = 0) . We evaluate this quantity 
in the following way UJJI. We use the double Mellin representation for the pion transition form factor reducing the 
integral in Eq. ([4]) to the convolution of propagator-like expressions. Then we perform the loop integration by using 
the standard Feynman a-representation. Finally, we are able to expand the integral over the ratios of the electron 



1 For completeness we give explicit expressions for the amplitude A (g 2 ) = A (g 2 ) — .A(O) for different regions of q 2 : ReA(q 2 ) = 
^[Li 2 (-y (ql)) + ^ + iln 2 (-y(g 2 ))], Im^(g 2 ) = 0, for? 2 < 0; and ReA(q 2 ) = -j^ da (-20), Im^(q 2 ) = 
~~jf^2) arct §[/% 2 )] i for < q 2 < 4m 2 . Here (3{q 2 ) = y/1 - Am 2 /q 2 , f3{q 2 ) = y/4m 2 /q 2 - 1, 6 = arctg[l//3(g 2 )] , and 
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Cl2(z) = — J dt In 1 2 sin(f/2)| is the Clausen's integral. 
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and pion masses to the characteristic scale of the pion form factor A by closing the Mellin contours in the appropriate 
manner and take the leading term of expansion. We arrive at the following representation: 
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where the constant xp (m) is defined by 
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with F^ytj* (t, t) being the physical pion transition form factor given in symmetric kinematics for space-like photon 
momenta t — Q 2 = —q 2 > 0. One has to note that the logarithmic dependence on the scale fj, appearing in Eq. (fT0|) 
as a result of the decomposition of the integral over the dimensional variable t into two parts is compensated by 
the scale dependence of the low-energy constant xp (/•*) displayed in Eq. (JTTJ) . The obtained representation defines 
the unknown subtraction constant in dispersion formula, Eq. ([7]) , via the pion transition form factor in a simple and 
transparent way. It is consistent with the result of the chiral perturbation theory [3, HE OS E3) Gil where xp (m) 
is the unknown low-energy constant. The result of Eq. (|10p also agrees with the conclusions made in 12J where the 
inequalities m e << << A were exploited. 

3. In order to estimate the integral in Eq. (JTTJ) , one needs to define the pion transition form factor in symmetric 
kinematics for space-like photon momenta. Since it is not known from the first principles, we will adapt the available 
experimental data to perform such estimates. Let us first use the fact that i 7 ! n - 7 « 7 » (t,t) < F^ 7 « 7 » (t, 0) for t > 
in order to obtain the lower bound of the integral in Eq. (111[) . For this purpose, we take the experimental results 
from the CELLO [lj| and CLEO [2(| Collaborations for the pion transition form factor in asymmetric kinematics for 
space-like photon momentum which is well parametrized by the monopole form 
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Note that Sq LEO «ml as predicted by simple vector meson dominance models (VMD) [1, H, [T2|, and is not far from 



the asymptotic prediction of operator product expansion (OPE) QCD [2l|: 
type of the form factor one finds from Eqs. (|10[) and (fTTjl that 
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-23.2 ±0.1. 
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Thus, for the branching ratio we are able to establish the important lower bound which considerably improves the 
unitary bound given by Eq. © 
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Now let us proceed further in this manner and assume that the monopole form is also a good parametrization for 
the form factor in symmetric kinematics 
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The scale S\ can be fixed from the relation for the slopes of the form factors in symmetric and asymmetric kinematics 
at low t , 
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that gives Si = sq/2. Note that similar reduction of the scale is also predicted by OPE QCD from the large 
momentum behavior of the form factors: sf PE — Sq PE /3 [2l|. Thus, the estimate for the amplitude in the limit 
q 2 — > can be obtained from Eq. (|13p by shifting the lower bound by a positive number which belongs to the interval 
[31n(2)/2,31n(3)/2]. One finds 



A( q * = 0)=- 3 -ln(^)-S= -21.9 ±0.3, 



(17) 
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that corresponds to the value \p ( m p) = O- 1 =t 0-3 of the low-energy constant xp taken at the scale fi — m p . Using 
the obtained prediction for the subtraction constant, one can evaluate the branching ratio 



Btn 



(6.23 ±0.09) • 10 



-s 



(18) 



This is 3 standard deviations lower than the KTeV result given by Eq. @. One can convert the experimental data 
to obtain the restriction on the scale parameter si in Eq. (|17p . Then for the amplitude at q 2 — estimated from the 
experimental data for Br^Y, (7r° — > e + e~) one finds 



A KTeV (q 2 



0) 



-18.6 ±0.9. 



(19) 



Since the amplitude in Eq. (|17|) depends logarithmically on the scale parameter si, one needs to reduce the value 
of si by a factor of larger than 4. However, it obviously contradicts the experimental data for the slope parameters 
discussed above. 

4. Let us compare our estimates with the results obtained in other approaches. We start with the QCD sum rule 
approach ,23|. There, the pion form factor of the transition process 7*7* — > ir° was found in the form 
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where v.c. are small corrections from the vacuum condensates and s^ CDsr is the so-called dual interval parameter 

0.7 GeV 2 in the original work [23J]. By using the relation for the form factor slopes in 
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one can identify the duality parameter with the CLEO parameter 



Then by using Eq. (TlTj) one finds 
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-21.7±0.1, 



(21) 
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that corresponds to the rescaling s QCDsi = ^s^ 00151 /e \ and is well suited to the interval in Eq. (flT]) . The corresponding 
branching ratio is shown in Table 1. 



TABLE I: Values of the quantity A (q 2 = 0) and the branching ratio B (tv° — > e + e ) obtained in our approach (see, Eq (|17[) ) 
and compared with various phenomenological models and the KTeV experimental result. 





CLEO+OPE 


QCDsr 


gVMD 


QM[12J 


N X QM 


NQM [10] 


Experiment [1] 




21.9 ± 0.3 


21.7 ±0.1 


21.9 


23.4 ±0.5 


22.1 ±0.5 


24.5 


18.6 ±0.9 


B (tt° -> e+e') x 10 s 


6.23 ±0.09 


6.21 ±0.05 


6.2 


5.8 ±0.2 


6.1 ±0.2 


5.38 


7.49 ± 0.38 



5. Then, we consider the parametrization of the pion form factor motivated by the generalized VMD [24 



-igVMD 



An 2 / 2 {s + t)st- h 2 st + h 5 (s + t) + M^M^h 7 
3 (M 2 + s) (M 2 + t) [Ml x + s) (M 2 X + t) 



(24) 



with the parameters M v = 769 MeV, M V1 = 1465 MeV, h 2 = -10 GeV 2 , h 5 = 6.93 GeV 4 , h 7 = 3/ (4tt 2 / 2 ). This 
parametrization satisfies above mentioned constraints on the pion transition form factor. At zero virtualities, the 
form factor is normalized by the axial anomaly. The relation in Eq. (|16[) is valid, which yields for the radius 



(r2)g VMD = _ g 



OF (t, 0) 



at 



= 6 



t=o 



M 2 



= 0.39 fm 2 



(25) 
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which is close to PDG average (r 2 )™^ = 0.407 ± 0.051 fm 2 0. Note also that numerically the second and third 

terms in Eq. (|25|) almost cancel each other. The form factor F*,J? (s, t) has also correct OPE QCD motivated 
behavior at large virtualities [2l| 



= ^i. (27) 

The constant /15 was fixed in [24[ from a fit of CLEO data (20j|, the constant /12 defining the next-to- leading order 
power correction to the form factor at large t was estimated in [25[ using the results of QCD sum rules [26j]. In the 
asymptotic OPE QCD limit, one has h 5 / (M V M V1 ) -> 6. 

For the generalized VMD form factor we estimate the subtraction constant 

*™ ( ,_ 0) _^ h (£) + . + ^_»*- hpf 



((r + 1) inr - 2 (r - 1)) - (l + ^) (, 2 - 1 - 2Hn,)) 



-21.94, (28) 



M v (r — l) 3 \2M V " V M V i M v. 

where r = (Af yi /Af y ) . It agrees well with our prediction given by Eq. (117|) . 

6. Let us now consider the amplitude 7r° — > e + e~ in the context of the constituent quark models. Within this 
model, the pion form factor is given by the quark-loop (triangle) diagram. Taking the constituent quark mass in the 
loop to be momentum independent, the result for the form factor in symmetric kinematics is given by [27| 



om 2M? / B„ (t) + 1 \ / Mj 

F?™, (t, t) = — i- In g , + , /3 g (t) = A/1 + 4-^. (29) 



q (t)t \P q (t)-l 

Substituting it into Eq. pT|) . one can get 



^QM(<Z 2 = 0)=31n^)-^, (30) 

in accordance with [jjj • For the constituent quark mass in the interval M q = 300 ± 50 MeV one finds 

A QM (q 2 = 0) --(23.4 ±0.5), (31) 

which is in agreement with the above estimates but contradicts the experimental result. In order to fit the KTeV value 
in Eq. (fl9|) . one needs to take the quark mass M q < 100 MeV, which is an unacceptable region for the constituent 
quark mass. 

By using the results for the pion transition form factor [28| obtained in the nonlocal quark model based on the 
instanton picture of QCD vacuum [2!|, one finds numerically from Eq. (|10[) that 



A xQM (q 2 = 0) = - (22.1 ±0.5), (32) 

which is consistent with Eq. (|17p . Within the instanton model, the quark mass is momentum dependent and for 
simplicity is taken in the Gaussian form 

M(fc) = M g exp(-A : 2 /A 2 QM ), (33) 

with M q = 300 ± 50 MeV and A x qm ~ 1 GeV fixed by fitting the pion decay constant f n . 

The similar result was obtained in in the nonlocal quark model (NQM) 0, ED] (see, Table 1). This model is 
based on the assumption that the quark propagator is described by the entire function without any singularities in 
the momentum space. Such an assumption guarantees the quark confinement in a sense that eliminates the threshold 
cuts corresponding to free quark production. A unified and uniformly accurate description of broad range of physical 
observables is obtained within the framework of this model, see, e.g. 31 1. 

One has to note that the nonlocal quark model with momentum dependent mass has an important advantage as 
compared to the mass independent (local) quark model. The first model has correct large momentum behavior given 
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by Eq. (|27|) . while the latter has an extra log(i) term in the asymptotic, as follows from Eq. (|29|) . For this reason, 
the local model is not suitable for fitting the CLEO data on the pion transition form factor. 

7. The 7] — > decay can be analyzed in a similar manner. As in the pion case, the CLEO Collaboration has 
parametrized the data for the 77-meson in the monopole form [20| : 



F^?(t,0)= } ChEQ , s^ = ( 774 ± 29 MeV ) Z > ( 34 ) 
1 -t- i/s 0v 

which is very close to the relevant pion parameter. Then following the pion case (with evident substitutions), one 
finds the bounds for the q 2 — > limit of the amplitude 77 — > n + fj," as 



o / CLEO\ r 

A W = 0) > -| ^ - j = - (7.2 ± 0.1) , (35) 



and for 77 — > e + e one gets again Eq. (|13p . The obtained estimates allow one to find the bounds for the branching 
ratios 

B{j] -> < (6-23 ±0.12) • 1CT 6 , (36) 

B (ji —f e+e") > (4.33 ± 0.02) • 10~ 9 . 

Note that for the decay 77 — > /i + /~t~ we get the upper limit for the branching. This is because the real part of the 
amplitude for this process taken at the physical point q 2 = m 2 for the parameter Sq ^ EO remains negative and a 
positive shift due to the change of the scale so n — > si v reduces the absolute value of the real part of the amplitude 
|Re^(m2)|. At the same time, considering the decays of tt° and 77 into an electron-positron pair, the evolution to 
physical point © makes the real part of the amplitude to be positive for the parameter Sq leo and the absolute value 
of the real part of the amplitude increases in changing the scales of the meson form factors. The predictions for the 
decays 7/ — > l + l~ obtained by reducing the scale s^ KO —> si v for the case of the 77-meson transition form factor are 
given in Table II. 



TABLE II: Values of the branchings B (P — » l + l ) obtained in our approach and compared with the available experimental 
results. 



B 


Unitary bound 


CLEO bound 


CLEO+OPE 


Experiment 


B (tt e + e~) x 10 8 


> 4.69 


> 5.85 ±0.03 


6.23 ±0.09 


7.49 ±0.38 [1] 


B (77 -> mV~) x 10 6 


> 4.36 


< 6.23 ±0.12 


5.11 ±0.20 


5.8 ±0.8 [7, 32J 


B (77 -> e+e") x 10 9 


> 1.78 


> 4.33 ±0.02 


4.60 ± 0.06 





8. In this work, we have derived in the leading order in (m e /A) the representation for the amplitude of the rare 
7T° — > e + e~ process in the limit q 2 — > 0. It is given in terms of the inverse moment of the transition pion form factor 
in symmetric kinematics for space-like photon momenta. By using data of the CELLO and CLEO Collaborations on 
the pion-photon transition form factor in the obtained representation, we found the new lower bound for the decay 
branching ratio which essentially improves the well-known unitary bound. Further constraints follow from the results 
of OPE QCD correlating the pion transition form factor in different kinematics as the change of characteristic scales. 
These considerations allow us to reconstruct the full decay amplitude and make predictions for the decay branching. 
A similar procedure is also applied to the decays 77 — > l + l~ . We compared our predictions with the results obtained 
in various phenomenological approaches and found that all of them are in agreement with our results. However, the 
obtained prediction for the branching ratio ir° — > e + e~ is 3a below the recent KTeV measurement. 
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